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Abstract 

In order to evaluate the response to cosmic-ray nuclei of a Compton-Scatter Transition Radiation Detector in 
the proposed ACCESS space-based mission, a hybrid Monte Carlo simulation using GEANT3 and an external 
transition radiation (TR) generator routine was constructed. This simulation was employed to study the effects 
of delta-ray production induced by high-energy nuclei and to maximize the ratio of TR to 5-ray background. The 
results demonstrate the ability of a Compton-Scatter Transition Radiation Detector to measure nuclei from boron 
to iron up to Lorentz factors 7 ~ 10 taking into account the steeply falling power-law cosmic ray spectra. 
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1. Introduction 

The proposed Advanced Cosmic-ray Composi- 
tion Experiment for Space Science (ACCESS) is 
a dedicated space-based mission to perform di- 
rect nuclear composition measurements up to near 
the 'knee' of the cosmic-ray spectrum, E ~ 10 6 
GeV/nucleus [1]. One configuration [2] employs a 
nearly 6 meter 3 Compton-Scatter Transition Ra- 
diation Detector (CSTRD) in a cubic geometry 
with a smaller calorimeter underneath. The abso- 
lute size of the ACCESS instrument is driven by 
the requirement to obtain significant event statis- 
tics in this flux-poor energy regime while being 
constrained by the achievable spacecraft size and 
mass. 



The Compton-scatter technique exploits the 
facts that 1) the x-ray spectra of emitted TR 
can be hardened by a judicious choice of radia- 
tor material, thickness, and spacing and 2) the 
greatest energy dependence of the TR is at high 
frequencies [2]. For a radiator with thickness l\, 
plasma frequency u>i, and spacing I2, interference 
effects from the superposition of amplitudes from 
each interface yield pronounced maxima and min- 
imum in the TR spectra. The highest frequency 

maximum is at uj m ax = — ^(1 + P) while the sat- 
uration energy is 7 S « 
P 

constants (i.e., metals) where p = 1. (In order to 
be conservative, the simulations employed in this 
study assume p = 0). Thus an ensemble of alu- 



^V l Ml + p), where 
except for materials with complex dielectric 
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minum foil radiators with 100 — 200 /im thickness 
separated by several mm leads to 7 S w 10 5 and 
fo^max ~ 130 keV (for p = 0). As the Compton- 
scattering cross-section is significant at these en- 
ergies, these photons will scatter away from the 
initiating particle's trajectory and thus spatially 
separate the TR signal from the ionization path. 
Furthermore, the significant flux of x-rays with 
energies > 50 keV can be detected with an in- 
organic (e.g., Csl) scintillator. The principle of 
Compton-scattcr TRDs has been demonstrated in 
accelerator test beams [3]. 

The thick radiators used to produce such hard 
x-rays also induce <5-ray production. The impact of 
this background is exacerbated in cosmic-ray mea- 
surements by the vast number of sub-TR events 
resulting from the power-law nature of the cosmic- 
ray spectrum. The study presented in this paper 
optimizes the CSTRD signal in relation to this 
background. 



2. Simulation Studies 

A hybrid Monte Carlo simulation has been 
constructed using GEANT3 [4] interfaced to an 
external TR generator routine based on the for- 
malism of Ter-Mikaelian [5]. GEANT3 models the 
dominant electromagnetic processes in the energy 
range 10 keV to 10 TeV with the cross-sections 
smoothly extending to higher energies. For heavy 
ions, ionization losses are simulated along with 
the discrete process of 5-ray production. Hadronic 
interactions are not simulated for heavy ions 
in GEANT3, and the ACCESS calorimeter and 
ancillary charge measurement detectors are not 
modeled in this study. 

Although multiple CSTRD configurations were 
investigated, the results employing a radiator ar- 
rangement with the highest TR saturation Lorentz 
factor are presented. The configuration consists of 
a large radiator stack surrounded on all six sides by 
highly segmented scintillators which measure both 
the Compton-scattered TR and ionization losses, 
thus providing the primary particle charge deter- 
mination. The radiator consists of 150 Al foils with 
150 fim thickness separated by 1 cm vacuum gaps 



Table 1. Response of a CSTRD to normally incident 
carbon nuclei with Lorentz factors above TR saturation. 
The first number in each column represents the number 
of pixel hits per event with 50 keV < E p i xe i < 500 keV 
while the second number (in parentheses) gives the pixel 
hits per event without energy selection. 

Csl layer TR only 5-rays only Sum TR +5-rays 

Top 6.7 (7.9) 0.5 (0.7) 7.1 (8.6) 7.1 (8.6) 

Sides 23.8 (27.9) 1.5 (2.3) 25.3 (30.2) 26.3 (31.4) 

Bottom 18.3 (22.3) 1.7 (4.5) 20.0 (26.8) 19.7 (26.3) 

Total 48.8 (58.1) 3.7 (7.6) 52.4 (65.7) 53.1 (66.3) 

with an areal size of 160 x 160 cm 2 . The radiator 
stack is surrounded by six 160 x 160 cm 2 planes of 
Csl scintillator in a cubic geometry. Each scintil- 
lator module contains a layer of 2 mm thick Csl, 
250 /im silicon PIN diode readout, and 5 mm of 
G10 encapsulated in a 1 mm Al housing. The Csl 
is segmented into 2.5 x 2.5 cm 2 pixels with each 
pixel recording the deposited energy on an event- 
by-event basis with a 50 keV threshold energy. 

The pixel TR energy distribution peaks at an 
energy per pixel near 100 keV with nearly 100% of 
the Compton-scattered TR signal being recorded 
below 500 keV. The <5-ray pixel energy distribution 
shows a low energy peak below 100 keV with a pro- 
nounced tail extending well above 1 MeV. A pixel 
energy selection of 50 < E pixe i < 500 keV retains 
virtually 100% of the TR hit pixels while rejecting 
> 40% of the 5-ray hits. The pixels recording the 
ionization of the primaries and their nearest neigh- 
bors are excluded in this pixel counting. 

Table 1 details the response for normally inci- 
dent carbon nuclei with 7 > 10 5 for simulated TR 
only, 5-ra.y background only, and signal + back- 
ground. Note that the pixel occupancy is at such a 
low level that summing the individual signal and 
background samples equals the result for simulta- 
neously simulating the two processes; i.e., there are 
very few events in which a <5-ray and TR photon 
hit the same pixel. The results indicate that the 
signal is 48.79/3.65 = 13.4 times larger than the 
background. By varying the scintillator thickness 
from 1 to 8 mm, the signal/background ratio was 
maximized at 2 mm Csl thickness. Furthermore, 
a substantial reduction in the 5-ray background 
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Fig. 1. The simulated total number of pixels hit for normally incident carbon (left) and iron (right) nuclei. The response 
to the background <5-ray process alone is illustrated by the lower, horizontal distribution. 



was obtained by placing the scintillators on the 
outside of the radiator stack versus interspersing 
them within the stack. No significant gain in the 
signal/background ratio was obtained by rejecting 
pixels further from the trajectory of the primary 
than just the hit pixel plus nearest neighbors. 

These simulations were performed using nuclear 
primaries from boron to iron and incident energies 
chosen from a flat log-E distribution to enhance 
the measured high energy response. Figure 1 shows 
the results for carbon and iron for TR + back- 
ground and for <5-ray background alone. This radi- 
ator configuration yielded a mean saturated value 
of - 960 TR photons for carbon and - 18, 000 TR 
photons for incident iron nuclei with, on average, 
~ 27% (20%) of the TR photons at energies greater 
than 50 (100) keV . 

The effective dynamic range of the TR mea- 
surement is obtained from a statistical analysis. 
Assuming a power law spectrum tp(E) ~ E~ a 
and a TR threshold energy of E t hr, the num- 
ber of events above the TR threshold goes as 
N signal ~ El^ r a . Some number of sub-threshold 
events accompanied by <5-rays will masquerade as 
high-energy events with TR. If the number of these 
misidentified low-energy events is required to be 



no more than f3N S i gna i, then the probability of a 
low-energy event with <5-rays fluctuating upward 
to mimic the signal of a valid event above E t hr 
must be no more than P = (3(Eq / 'Ethr) ' 1 ■ Using 
the conservative values of a = 2.75, (3 — 0.01, 
E = 1 GeV, and E thr = 10 5 GeV, the probability 
is evaluated as approximately 2 x 10~ n or ~ 7 a 
assuming a single-sided Gaussian probability dis- 
tribution function. Thus the low end of the TR 
dynamic range is defined as that energy where 
the TR response is 7 a above the average value of 
the sub-TR response. On a power-law decreasing 
spectrum, the flux of particles at the high end of 
the energy range is depleted as compared to lower 
energies. Thus, the high end of the TR dynamic 
range can have a more modest 2 a separation 
between the signal and the saturated TR value. 

Table 2 details the effective dynamic ranges for 
various nuclear species based upon the results of 
this analysis. As the data demonstrate, the mod- 
eled CSTRD can measure the Lorentz factors of in- 
cident nuclei from boron to iron. Furthemore, the 
range of Lorentz factors that guarantee signal sep- 
aration from the sub-TR background and the sat- 
urated value significantly increases as the atomic 
number of the incident nuclei increases. 
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Tabic 2. Dynamic range of TR measurements for 
various incident nuclei. 



Nuclei Effective n ow Effective ■Jhigh Dynamic Range 



11 Boron 


2.4 X 10 4 


4.5 x 10 4 


1.9 


12 Carbon 


2.2 X 10 4 


5.2 x 10 4 


2.4 


16 Oxygen 


1.6 x 10 4 


5.8 x 10 4 


3.6 


28 Silicon 


1.1 x 10 4 


7.2 x 10 4 


6.5 


56 Iron 


7 x 10 3 


8.3 x 10 4 


11.9 



Although hadronic interactions of the incident 
nuclei are not modeled for this study, an estimate 
of the attenuating effects of the relatively thick ra- 
diators and Csl scintillators in the CSTRD can be 
obtained by consideration of the nuclear inelastic 
cross-sections. For iron primaries incident upon an 
aluminum target, the nuclear interaction length is 
given as 20.3 g/cm 2 [6] while the grammage of the 
CSTRD modeled in this study is 11.5 g/cm 2 . Thus, 
approximately 43% of an incident iron flux will 
be lost due to hadronic inelastic collisions. Events 
with charge-changing interactions will be recog- 
nized by observing a difference in the charge mea- 
surement of the particles entering versus exiting 
the detector. 
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3. Conclusions 

A Monte Carlo simulation has been employed to 
simulate the response of a Compton-scatter tran- 
sition radiation detector and to optimize the de- 
sign with respect to the background process of 5- 
ray production. The results indicate that a space- 
based CSTRD can perform nuclear composition 
measurements up to a Lorentz factor of 10 5 . How- 
ever, while the dynamic range is more than a factor 
of 10 for cosmic-ray iron measurements, the range 
is more limited for lower- Z nuclei due to the rela- 
tive reduction in TR generation. 
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